Effects of thermal fluctuations on the magnetic behavior of mesoscopic 

superconductors 
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We study the influence of thermal fluctuations on the magnetic behavior of square mesoscopic 
superconductors. The strength of thermal fluctuations are parameterized using the Ginzburg num- 
ber, which is small (G, ~ 10 -10 ) in low-T c superconductors and large in high-T c superconductors 
(Gi » 10~ 4 ). For low-T c mesoscopic superconductors we found that the meta-stable states due to 
the surface barrier have a large half-life time, which leads to the hysteresis in the magnetization 
curves as observed experimentally. A very different behavior appears for high-T c mesoscopic su- 
perconductors where thermally activated vortex entrance/exit through surface barriers is frequent. 
This leads to a reduction of the magnetization and a non-integer average number of flux quanta pen- 
etrating the superconductor. The magnetic field dependence of the probability for the occurrence 
of the different vortex states and the fluctuations in the number of vortices are studied. 

PACS numbers: 74.78.Na, 74.40.+k, 74.20.De 



I. INTRODUCTION 

In the last decade, the study of vortex matter 
in finite mesoscopic superconductors attracted a lot 
of attention, both experimentally 1 : 2 : 3 : 4 : 5 : 6 : 7 ! 8 ! 9 ! 10 and 
theoretically. 11 ! 12 ! 13 ! 14 ! 15 : 16 ! 17 ! 18 ! 19 Most of the experi- 
mental studies on mesoscopic superconductors deal with 
conventional low-T c superconductors. Resistivity mea- 
surements are widely used to obtain information on the 
superconducting/normal transitioniS^l With ballistic 
Hall magnetometry it is possible to measure deep inside 
the superconducting phase diagram and to obtain infor- 
mation on the magnetization of the different supercon- 
ducting statesi 5 ! 6 : 7 : 8 Recently, A. Kanda et al. devel- 
oped the multiple-small-tunncl-junction method to dis- 
tinguish between multivortex and giant vortex states 
in mesoscopic superconductors Aifi Theoretically, stud- 
ies of the Ginzburg-Landau free e nergjiii4i2ii 3 -i 4 *i 5 jiSii 7 . 
and the time dependent Ginzburg-Landau equations 1 ^ 9 - 
have shown how the magnetic and dynamics properties 
depend on the sample sizes and geometry. In particular 
the mesoscopic samples can develop Abrikosov multivor- 
tex states^ and depending on the size of the sample it 
is possible to observe first or second order transitions. 11 
One interesting characteristic of the magnetic proper- 
ties of mesoscopic superconductors is the behavior of the 
dc magnetization curves. In mesoscopic superconductors 
there is a reinforcement of the surface barrier for entrance 
and exit of vorticesii 3 ^* 1 ^ The surface barriers allow for 
the existence of meta-stable states of constant vorticity 
as a function of magnetic field and lead to a magneti- 
zation curve with discontinuous jumps, which have been 
observed experimentally^ All these results are for con- 
ventional low-T c superconductors, where the effects of 
thermal fluctuations are generally small. 

A very different situation occurs in macroscopic high- 
Tc superconductors where the effects of thermal fluc- 



tuations are important. In particular, fluctuations af- 
fect the dynamics of vortex entrance through the sur- 
face barriers i^SiSiiSSiS^iSii Recently, there has been interest 
in the study of the magnetic properties of micron-sized 
high-r c superconductors^ 5 ^. Single fluxoid transitions 
and two-state telegraph noise was measured in thin film 
rings at temperatures close to T c m& The fluxoid transi- 
tions were explained by thermally activated penetration 
of vortices. 

The magnetic and dynamic behavior of mesoscopic su- 
perconductors is greatly influenced by surface barriers 
and geometry effects ^iiiiLi^ As a consequence, ther- 
mal fluctuations in mesoscopic high-T c superconductors 
could lead to interesting phenomena due to thermally 
activated vortex entrance/exit through surface barriers. 

In the present paper we study the effects due to ther- 
mal fluctuations on the vortex entry in square meso- 
scopic superconductors. We solve numerically the time- 
dependent Ginzburg-Landau equations taking into ac- 
count demagnetizing effects and thermal noise fluctua- 
tions. The strength of thermal noise fluctuations is pa- 
rameterized by the Ginzburg number. For conventional 
low-Tc superconductors the Ginzburg number is small, 
while for high-T c superconductors this parameter is large. 
We compare results obtained from the minimization of 
the mean field free energy, corresponding to the equilib- 
rium states in the absence of thermal fluctuations, with 
the results from the time dependent Ginzburg Landau 
(TDGL) dynamics with small Ginzburg number (corre- 
sponding to conventional low-T c superconductors) and 
the results for TDGL dynamics with large Ginzburg num- 
ber (corresponding to high-T c superconductors). The 
magnetization, the number of vortices, the free energy, 
etc are calculated in each case. We find very good agree- 
ment between the results of free energy minimization and 
the results of TDGL with small Ginzburg number. In the 
case of a large Ginzburg number the TDGL approach 
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gives different and new results for the magnetization and 
the vortex dynamics. 

The paper is organized as follows. In Section II we 
give the two theoretical formalisms used in the present 
paper. First we describe how we solve the TDGL equa- 
tions, obtaining the magnetic fields via the Biot-Savart 
law and including thermal fluctuations. Secondly, we de- 
scribe the solution of the stationary Ginzburg-Landau 
equations, when the magnetic field is solved via Ampere's 
law. In Section III the results of the TDGL theory are 
compared with the ones of the stationary GL theory. In 
a first step, we neglect the thermal fluctuations and com- 
pare the results both at zero and finite temperature. In 
a second step, we include thermal fluctuations in our cal- 
culations and compare the results with small and large 
Ginzburg number. Finally, in Section IV we present our 
conclusions. 



II. THEORETICAL FORMALISM 

A. Time-dependent Ginzburg-Landau theory 

We solve the time-dependent Ginzburg-Landau equa- 
tions, in the gauge where the electrostatic potential is 
zero, and t aking into account thermal fluctuations (see 
Ref. 127128ft . i.e. 
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where VP is the order parameter and A the vector poten- 
tial, r^, = 2mD/h 2 and — c 2 /<r n are kinetic coeffi- 
cients where a n is the quasiparticle conductivity and D is 
the electron diffusion constant. Q<s[f, t) and ^(r, i) are 
Langevin thermal noise with average (f$r) = ((a) = 
and the following correlations 2 ^* 2 ^: 

(a(f f ,t)^(f,t'))=2k B Tr^d(f-f^)d(t-t'), (3) 
(GfatXi^O) = 2k B TT?8(f- ?)6(t - (4) 

The difference Q between the superconducting and the 
normal state Gibbs free energy is given by 



G = Gs — Gn 

a(T)|#| 2 



2m, 



-L(VxA-H a ) 
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d 6 f. 



(5) 



To include demagnetization fields, the magnetic term in 
G (last term) must be integrated inside and outside the 
sample volume. The differential Eqs. (JTJ and J2J) must 
be solved in the whole space. 



Taking the variations in Eqs. (1) and (2), the normal- 
ized TDGL equations becomo 18 i 19 i 27 i 29 i 30 
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Lengths have been scaled in units of the coherence length 
£(0), times in units of to = Aira n \ 2 L / 'c 2 , A in units of 
# c2 (0)£(0), * in units of #«, = [mc 2 /87re 2 A(T) 2 ] 1 / 2 and 
temperature in units of T c . rj is equal to the ratio of the 
characteristic time to for the relaxation of A and the time 
tcL for the relaxation of \&: rj = tcL/to = c 2 /(47rcr n K 2 D), 
with toL = £ 2 /D. For superconductors with magnetic 
impurities we have Di mp = c 2 /(487T(t„k 2 ), and therefore 
i] = 12 in this case. 

The normalized thermal noise correlations are 

(& (f, *)C*(**.0) = —V2GlT5(f-r<)5(t-t'), (8) 
V 

(CW, t')) = A-K^2G{T5{r- ?)8(t - (9) 

where 



Gi = 



1 



H 2 (o)e(o) 



(10) 



is the Ginzburg number- 1,32 that governs the strength of 
thermal fluctuations. Gi measures the relative size of the 
minimal (T=0) condensation energy H 2 (0)£ 3 (0) within a 
coherence volume and the energy of thermal fluctuations 
at T c , ksTc. For conventional low-T c superconductors 
Gi « 10~ 8 and for high-T c superconductors Gi w 10~ 4 
(see Ref. HI). 
Equation J7J is Ampere's law, i.e. 



J = k 2 V x V x A, 
where J is the current inside the sample: 
J(r) 

J(r) 



= J n (r) + J s (r) + ( A (r,t), 
_dl 
~~dt 
+ U(r,t). 



(l-T)Im[**(V-M)*] 



(11) 
(12) 

(13) 



J n is the current due to the normal electrons, J n 
and J s is the supercurrent. 

The Biot-Savart law relates the magnetic induction 
B = V x A, the applied magnetic field H a and the sample 
currents: 



B(r) - H a 
B(f) - H a 



1 



Attn 2 



J(r') x 



r'|3 



Q(r,r')g(r')d S r'. 



(14) 
(15) 
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This equation states that the magnetic field is completely 
determined by the externally applied field H a and by the 
currents flowing inside the sample (which are given by 
g{r')). The scalar function g(r') is the local magnetiza- 
tion or density of tiny current loops and is given by2£ 



J{x,y) = -z x Vg(x,y) = V x zg(x,y). 



(16) 



This guarantees that V • J = 0. 

The kernel Q can be obtained from the Biot-Savart 
law. In the dipolar approximation, i.e. for long distances, 
p = \f- r'\ and Q = -l/Airp 3 . 

Then g(r) and Q(r, r') have all the information that we 
need to include the demagnetization fields (and Q(r, r') 
is a known and time- independent kernel). If we write 
the TDGL equations in terms of g(f), and we use the 
Biot-Savart law, we only need to solve the differential 
equations inside the sample volume and the demagneti- 
zation fields are automatically included^ On the other 
hand, we can also solve Ampere's law retaining the 3D 
nature of the magnetic field distribution. Both meth- 
ods are equivalent for steady-state magnetic phenomena 
(V • J = 0) when the Biot-Savart and Ampere's laws are 
equivalent^ 

In the present paper we will deal with square samples 
with thickness d (d < A). In this case, the order param- 
eter and g(x,y) can be assumed to be uniform in the z 
direction. Using the Biot-Savart law it is only necessary 
to write the TDGL equations for a 2D sample^ i.e. 



~dt 



dA 
dt 



-[{I - T)*(l - |*| 2 ) - {-iV 2D - I) 2 *] 



(17) 



B z (x,y) 



(1 - T)Im[**(V 2 _D - iA)% - V x zg 

a (is) 

H a = [ Q(R,R')g(R')d 2 R', (19) 

J A 



were R = (x,y) and R' = (x',y'). We assumed H a \\z. 
To obtain g(x,y) we invert Equation l|19|) at each time 
step using the efficient Conjugate Gradient Method^ 

The boundary conditions for this problem, at the sur- 
face of the superconducting sample, are 



(-iV 



2D 



A)*\ n = 0, 
Sib = o. 



(20) 
(21) 



The system of GL equations, using dimensionless vari- 
ables and the London gauge V • A = 0, has the following 

form> ll i 12 i 17 



(-iVao -£) * = * (l - |*| 



-A 3D A = — S(z)j 2 i 



(22) 
(23) 



where 



32D = — (**V 2£ >* 



fV 2 D$*)-|*l A ( 24 ) 



is the density of superconducting current. The order pa- 
rameter satisfies the boundary conditions: 



-iV 



2D 



A) * 



= 0. 



(25) 



The boundary condition for the vector potential in carte- 
sian coordinates becomes: 



A\\ x \=R a ,\ y \=R, 



H {x,-y)/2, 



(26) 



at the boundary R s of a larger space grid. 

Here the distance is measured in units of the coherence 
length ^, the vector potential in cfi,/2e£, and the magnetic 
field in H c2 = ch/2e£ 2 = Ky/2H C . The superconductor is 
placed in the (x, y) plane, the external magnetic field is 
directed along the z axis, and the indices 2D, 3D refer to 
two- and three-dimensional operators, respectively. 



III. RESULTS 

We will compare the results obtained from minimiza- 
tion of the GL free energy (with magnetic fields solved 
via Ampere's law) with the stationary states obtained 
from the dynamics of the TDGL equation (with mag- 
netic fields solved via Biot-Savart law). We study the 
influence of the square symmetry restrictions on the vor- 
tex entrance in square samples. In a first step we will 
neglect the thermal fluctuations and we will investigate 
vortex states in mesoscopic superconducting squares at 
T = and T = 0.7T C , both within the time-dependent 
GL theory and within the stationary GL theory. In a 
second step we will take into account the thermal fluc- 
tuations in the TDGL theory and again we will compare 
the obtained results with the ones from the stationary 
GL theory. Such a comparison will be performed both 
for conventional low-T c and for high-T c superconductors. 
In the last part we focus on the time average of the order 
parameter l 1 ?!^ and the number of vortices N v for low 
and high-T c superconductors. 



B. Stationary Ginzburg-Landau theory 

To calculate the equilibrium states we minimize the 
mean-field free energy without thermal fluctuations and 
retain the three-dimensional magnetic field distribution. 



A. No thermal fluctuations 

First we consider superconducting squares with Gi = 
0, for which, thermal fluctuations are not taken into ac- 
count. Fig. n snows our results as a function of the ap- 
plied magnetic field for a superconducting square with 
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sides equal to W = 10£(0) and thickness d — 0.1£(0). 
The GL parameter equals k — 1 and temperature T = 0. 

Throughout the paper, we will show results of the ap- 
parent magnetization defined as 47rM a = (B z ) — H a . The 
real magnetization is equal to AttM = (B z ) — H where 
H is the internal magnetic field that can also be written 
as H = H a — (AttM)N where N is the demagnetization 
factor that depends on the sample geometry. However, 
M a is what is measured experimentally. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 



H a /H c 2 ( T ) 

FIG. 1: Results at G, = and T = 0. (a) The squares 
(■) show the apparent magnetization corresponding to the 
energy ground state and the straight lines the M a results of 
all the meta-stable states, (b) Gibbs free energy of the differ- 
ent vortex states, the inset is a zoom at low fields. Apparent 
magnetization (c) and the number of vortices (d) calculated 
with the TDGL equations, when increasing (a) and decreas- 
ing (V) the external magnetic field. The arrows show the 
vortex transition with TDGL. 

The square symbols in Fig.QIa) give the results of the 
apparent magnetization, corresponding to the "ground 
state" (global minimum of the mean-field free energy) as 
calculated within the stationary GL theory. The ground 
state vortex transitions are found calculating the free en- 
ergy of all the possible (meta-stable) vortex states that 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



H /H ,(T) 

FIG. 2: Similar results as Fig. □ but now at T = 0.7T C 
and Gi — 0. (a) Apparent magnetization, M a . The results 
obtained minimizing the mean-field GL free energy (■) are 
compared to the dynamical results calculated with the TDGL 
equations, when increasing (a) and decreasing (V) the exter- 
nal magnetic field, (b) The total vorticity vs. H a for the same 
cases as in (a), (c) Gibbs free energy of the different vortex 
states as a function of H a - 

are shown in Fig. ^b). The inset is a zoom of the 
free energy at low fields. The different branches cor- 
respond to different values of the total vorticity, N v = 
§{A + T^j5)rfi. Then, according to the global minima 
of the free energy, the vorticity changes always by one 
unit (AJV„ = 1). 

In Fig.^c) we study the behavior of the apparent mag- 
netization within the framework of the TDGL theory. 
The triangles show the behavior for increasing (A) and 
decreasing (v) field. Notice that the transitions while 
increasing the field or decreasing the field do not occur 
at the same field values, i.e. the penetration field (H p ) 
differs from the expulsion field (H e ). This hysteresis is 
due to the surface barrier With increasing and de- 
creasing field we see that the vorticity changes by four 
and minus four at the transitions, which means that the 
vortices enter or leave the sample by four, i.e. one vor- 
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tex through each side of the square. This is shown more 
explicitly in Fig. dd) where the total vorticity, i. e. the 
number of vortices inside the sample, are plotted as a 
function of the applied magnetic field. Of course, follow- 
ing the ground state free energy, the vorticity changes 
always by one unit. The fact that in the TDGL results 
the vorticity changes by four at the transitions means 
that, in this case, there is not a simple dynamical path 
to break the square symmetry in the system during the 
field sweep. The close circles (•) in Fig. ^c) are the 
—M a / H c2 (T) values corresponding to the meta-stable 
states with vorticity equal to a multiple of four as cal- 
culated within the stationary GL theory. Notice that the 
values of —M a /H C 2(T) obtained using the TDGL and 
stationary GL theory are similar for the same number of 
vortices. 

Next, we investigate the influence of finite tempera- 
ture neglecting the thermal fluctuations, i.e. G,; = 0. 
Figs. EJa-c) show the same as Figs, ^b-d), but now at 
temperature T = 0.7T C . Whe n we incr ease the tem- 
perature, both £(T) = £(0)/y/l -T/T c and X(T) = 
A(0)/-\/l — T/T c increase and the size of the sample de- 
creases measured in units of £(T), (£(T = 0.7T c )/£(0) = 
1.83). In this case less vortex penetration events are nec- 
essary to arrive to the normal state. Notice that this is 
what happens in a typical experimental situation when 
the temperature is increased. The main difference with 
the zero-temperature situation is that in the present case 
the vorticity changes with AN V = 2 for increasing mag- 
netic field while the vorticity changes with A7V„ = — 1 
in most transitions for decreasing H a - This behavior is 
different from the results obtained in Fig. 1(d) because 
of the small size of the sample as compared to £ (T) . For 
example, in Fig. EJb) we plot the number of vortices for 
increasing (a) and decreasing ( V) the magnetic field and 
we compare these results with the vorticity of the energy 
ground state (■). The arrows show the vortex transitions 
with TDGL. In Fig- E^b) we see the presence of hysteresis 
due to meta-stable states. 



B. Thermal fluctuations 

Next, we include thermal fluctuations, i.e. G; =/= in 
Eq. GJJ. 

It is known that the free energy minima results ob- 
tained at temperature T\ can be rescaled and will be 
equal to the free energy minima results of an "equiva- 
lent" system at a different temperature T 2 ^ T%. To be 
equivalent the systems must have the same size measured 
in units of £(T) and the same k. In Fig. OJa) we com- 
pare the ground state apparent magnetization results for 
two situations: (i) a square with sides W = 10£(0) and 
thickness d — 0.1£(0) at T — 0, and (ii) for a square 
with W = 10£(T) and d = 0.1£(T) at T = 0.7T C for 
Gi = 1CP 10 . As the size of the system is the same mea- 
sured in £(T) the latter results can be compared if we nor- 
malize the field by the temperature dependent H c2 {T), 




-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
H/H ? (T) 



FIG. 3: Results obtained including thermal fluctuations for 
T = 0.7T C . (a) G z = 1(T 10 , (b) G % = 1CT 8 and (c) = 1CT 5 . 
The system size in (a) is W = 10£(T), d = 0.1£(T) and in 
(b) and (c) W = 10£(0), d = 0.1£(0). The triangles are for 
increasing (a) and decreasing (V) the magnetic field and the 
solid squares are the results corresponding to the mean field 
ground states. 



as we did in Fig OJa). We see that thermal fluctua- 
tions have broken the symmetry restrictions for vortex 
entrance and in this case A7Y„ = 2 for H a < 0.8H C 2(T) 
and AN V = 1 for H a > 0.8H c2 (T), while AN V = 4 when 
neglecting the thermal fluctuations (see Fig. ^a)). Even 
when AiV^ = 1, the dynamic TDGL results for — M a 
(open symbols) do not follow in detail the — M a curves 
obtained from free energy minimization (solid symbols). 
This behavior will be more clear in Fig.^b) in case of a 
smaller sample. 

In Fig. Gib) we show the apparent magnetization as a 
function of the applied magnetic field for a square with 
sides equal to W — 10£(0) and thickness d = 0.1£(0) at 
T = 0.7T C for d = 10~ 8 , i.e. a conventional low-T c su- 
perconductor. The square symbols and the solid curves 
are the apparent magnetization for the ground state and 
the meta-stable states as calculated within the framework 
of the stationary Ginzburg-Landau theory. The trian- 
gles indicate the results with increasing (a) and decreas- 
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ing (v) field when using the time-dependent Ginzburg- 
Landau theory, taking into account the thermal fluctua- 
tions. For this system size and Gi = 10 -8 the TDGL re- 
sults with increasing and decreasing magnetic field show 
that AN V = 1. But, even though AN V — 1, there is still 
hysteresis in the — (M a ) curve for such a conventional 
mesoscopic low-T c superconductor. The system clearly 
does not follow the mean field free energy minima be- 
havior. Notice that there is a rather good agreement 
between the results from the TDGL theory and the sta- 
tionary GL theory, both for the values of — M a and for 
the transition fields. 

Even when the thermal fluctuations are weak, as in 
Fig-EIb), it is expected that for a sufficient long simu- 
lation time the apparent magnetization would relax un- 
til the ground state is reached. In this case, the relax- 
ation of the magnetic flux is related with the possibil- 
ity that vortices overcome the surface barrier by ther- 
mal activation. In our simulations we allow the system 
to relax during 250000 time steps with At = 0.0045i o , 
which gives the system the possibility to relax during t = 
1125 to ~ 0.1/isec. However, the hysteresis in the appar- 
ent magnetization curve is also observed experimentally 5 
in mesoscopic low-T c superconductors. Consequently, we 
can conclude that the half-life time of the meta-stable 
states is very large in low-T c superconductors even in 
comparison with the usual experimental times. 

Next, we repeat the calculation for the same sample 
as in Fig- E^b) , but now including thermal fluctuations 
of strength corresponding to a high-T c superconducting 
material (Gi = 10 -5 ). The results are shown in Fig. Etc). 
We see that the first two penetration fields obtained from 
the TDGL equations agree very well with the ones re- 
sulting from the mean-field free energy minima. This 
means that the effect of the surface barrier is almost 
suppressed here, once we allow the system to relax a 
short time. On the other hand, the strong thermal fluc- 
tuations decrease the —(M a ) values in Fig. Etc), with 
respect to the mean-field values. Below the first pen- 
etration field, in the Meissner state, the magnetic flux 
inside the sample (BW 2 ) depends only on the penetra- 
tion length. Therefore, in the Meissner branch of figure 
Fig. Etc) we can deduce from the slope of M a (H a ) what 
appears to be an increase in A induced by thermal fluc- 
tuations. In fact, it is known that thermal fluctuations 
can induce an increase of the effective London penetra- 
tion length in high-T c superconductors (see for example 
which is a consequence of the renormalization of 
the superfluid density due to thermal fluctuations, since 
A| ff = mc 2 /(87re 2 <|#| 2 )). 

In order to understand these effects in Fig. Etc), we in- 
vestigate the time average of the order parameter (l^l 2 ) 
and the time variations in the number of vortices (N v (t)) 
for low and high-T c superconductors, with Gi = 10~ 8 
and Gi = 10 -5 , respectively. In Figs.Et a ) and^Jc) we 
show (l^l 2 ) along the x direction, i.e. from the mid- 
dle of one side to the middle of the opposite side of 
the square, for two values of the applied magnetic field, 



1.0 




2 4 6 8 10 -3 -2 -1 01234 



X/ X(0) 



0.8 




XU(0) N v 

FIG. 4: Time average of the order parameter in the x di- 
rection, |*|^, for d = 10" 8 (■) and d = 10" 5 (□) for (a) 
H a = 0A2H c2 (T) and (c) H a = 0.75ff c2 (T). (b) and (d) 
are histograms of the different vortex states obtained when 
G t = HT 8 . 



H a = 0A2H c2 (T) and H a = 0.75H c2 (T). We obtain 
the following results: for Gi = 10~ 8 the number of vor- 
tices do not change in time and the (l^j )(x) results (■) 
are near to the mean field ones (Gi = 0). However, for 
d = 10~ 5 we observe that (|*| 2 )(a;) (□) is lower than 
the mean field (MF) result and the number of vortices 
fluctuates in time. 

The field H a = 0A2H c2 (T), corresponds to the Meiss- 
ner state without vortices in the mean-field case, see 
Fig.Etb). For G l = 10~ 8 we also obtain that N v (t) = 
for all t. However, for Gi = 10 -5 the number of vor- 
tices change in time. In Fig. 0tb) we plot a histogram of 
the different vortex states obtained during 250000 time 
steps. We see that the state with N v = is the most 
probable one, but that there is also some probability to 
have jumps to other vortex states, mainly with N v = ±1. 
Thus the reduction of (| V I / | 2 ) from the mean field result 
for Gi — 10~ 5 is a consequence of thermal fluctuations 
which allow the nucleation of thermally induced vortex- 
antivortex pair o 38i39 ' 4 ° i41 as well as the entrance and exit 
of thermally activated vortices. This leads to an increase 
of the effective A, due to the lowering of (| V I'| 2 ). 

Results similar to Figs. Et a )-(b) are obtained in 
Figs. Etc)- (d) for H a = 0.75H c2 (T). Fig-Hc) shows that, 
for Gi = 10 -8 (■), we have a vortex state with one vor- 
tex: (l^l 2 ) = in the center of the sample where the 
vortex is located. However, for Gi = 10~ 5 (n), we ob- 
serve that (|^I'| 2 )(a;) is different from zero in the center 
of the sample. The corresponding vortex state histogram 
is shown in Fig. Etd). We see that the state with one 
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FIG. 5: (a) Time average of the different vortex states ob- 
tained for G l = 10" 5 for increasing (a) and decreasing (V) the 
magnetic field. The black squares (■) are the ground state 
results, (b) Probability of the different vortex states. The 
inset is a zoom at low fields in a logarithmic scale. 

vortex is the most probable one but that there is also a 
finite probability for other vortex states which contribute 
differently to (\^\ 2 ). This makes (\V\ 2 ) ^ in the center 
of the sample. There are continuous entrances and exits 
of vortices and the time average is a mixture of states 
with different vorticity. 

More details about thermal activation of vortices can 
be obtained by calculating the average number of vortices 
(N v ): 

At tf 

(N v ) = - r (27) 

tf u t=u 

and the probability to have N v = n vortices (P(Ar„ =n )): 

At tf 

P(N v =n) = V] $N v (t);n > (28) 

f 1 t=n 

both magnitudes can be calculated if we know the num- 
ber of vortices in each time step N v (t) . Where the sum 
in Eqs. (|27l) and i)28|) are taken between two time values 
(ti and tf), and 5^ (t);n is the Kronecker delta. 

In Fig. [5f a) we plot the average number of vortices 
and in Fig.[5{b) the probabilities to find a specific vortex 
state as a function of the applied magnetic field. These 
results where obtained during the magnetic field sweep 
up shown in Fig. Etc) (G, = 1CT 5 ). 



FIG. 6: (a) Fluctuations in the number of vortices and (b) 
in —M a / H C 2(0) for Gi = 10 -5 corresponding to the situation 
of Fig. 



In Fig. |S{a), for Gi = 10 5 (A) and small mag- 
netic fields {H a < 0.7H c2 (T)), we see that (N v ) vs. 
H a follows the mean field result [N^ F (H a )) (■). For 
H a > Q.7H C 2(T) the results are different, except at the 
start of each step where (N v ) = Nff . For strong ther- 
mal fluctuations we see that the average number of vor- 
tices can change continuously and is frequently larger 
than7V„ MF . 

More details about the thermal activation of vortices 
can be obtained calculating the probability to have N v 
vortices using Eq. (|28|l . In Fig. OJb) we show the proba- 
bilities of the different vortex states for Gi = 10~ 5 . The 
inset is a zoom of the low magnetic field region in a log- 
arithmic scale. Near H a = 0, the probability of the zero 
vortex state is near one (P(n v =o) ~ 1) an d we also see 
that P(x v —u pa P(n v =-i) ~ 0.01. From the inset we find 
that for H a w 0.5H C 2(T) there is a discontinuous jump 
down in P(jv l;= o) and that at the same time the proba- 
bility of having N v = 1 increases abruptly. When this 
happens we see that (N v ) = 1 in Fig.[3fa). 

We also calculated the fluctuations in the vorticity 
(N v (t)) and the fluctuations in -M a (t)/ H c2 (0). The fluc- 
tuations in the magnitude of the apparent magnetization 
M a (t) is obtained through the relation: 

AM a = ^(Ml) - (M a )* , (29) 

where the mean values of M a (i.e. (M a ) and (M 2 )) are 
calculated using the definition given in Eq. I|27|) . 

The fluctuations are shown in Fig. |SJ From Fig. EJa) 
we see that the fluctuations in the vorticity increase with 
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increasing field until the first penetration field, H p i , (first 
dashed vertical line), where there is a remarkable de- 
crease in the number of vortex fluctuations. The fluc- 
tuations in M a also decrease at H p \ as we can see from 
Fig-lUb). At the second penetration field (second dashed 
line) in Fig. H3 the decrease in the fluctuations at H P 2 is 
almost completely buried under the thermal noise due to 
a continuous increase in the fluctuations of the number 
of vortices with increasing field, as seen in Fig. IBTa). 

IV. CONCLUSIONS 

We studied the influence of the strength of thermal 
fluctuations on the magnetic behavior of mesoscopic su- 
perconductors. In the absence of thermal fluctuations 
and surface imperfections, we found that the vorticity 
changes by four at the transition fields due to the sym- 
metry of the mesoscopic square sample. The surface bar- 
riers and the symmetry of the sample produce hysteresis 
and meta-stable states. 

In low-T c mesoscopic superconductors, where thermal 
fluctuations are small, we did not find thermally acti- 
vated entrance/exit of vortices through surface barriers. 
This result agrees with experiments in low-T c supercon- 
ductors where hysteresis in the penetration fields and 
meta-stable states were found. 

In low-T c superconductors, small thermal fluctuations 
however are enough to break the square symmetry re- 
striction to vortex entrance and less vortices penetrate 
at the same field. 

A different behavior was found when the strength of 
thermal fluctuations is increased as is the case for high-T c 
superconductors. In this case there are frequent ther- 
mally activated events of entrance/exit of vortices in 
agreement with recent experimental results on micron- 
sized high-T c superconducting thin film rings 26 . We 
also found that thermal fluctuations increase the effec- 
tive London penetration depth if we compare it with the 
mean field result. 

In high temperature superconductors, the d-wave sym- 



metry of the ground state and the tetragonal symme- 
try of the underlying crystalline structure can have an 
important effect in some of their macroscopic proper- 
ties. These have been modeled by phenomenological 
Ginzburg-Landau or London theories, containing mixed 
gradient couplings to an order parameter with a differ- 
ent symmetry 4 ^^ or additional quartic derivative non- 
local terms^^ An important consequence of this is 
that, at large magnetic fields, the equilibrium vortex lat- 
tice deviates from the triangular Abrikosov lattice to a 
square lattice, as it has been observed recently^ Also 
the possible effects in the individual structure of vor- 
tices in cuprates superconductors has been the subject 
of intense research^ 7 . In the case of mesoscopic high-T c 
superconductors, the structure of the vortex configura- 
tions confined within the geometry of the sample and the 
quantitative values of the surface barrier and penetration 
fields could be affected by such considerations. On the 
other hand, the qualitative physics of the thermally ac- 
tivated entrance/exit of vortices through surface barriers 
should remain similar to the one described here with the 
dynamics of the simple GL equations. In any case, more 
experimental and theoretical works in this field are nec- 
essary. Therefore, the study of the physical properties 
of mesoscopic high-T c superconductors can contribute to 
the global understanding of these interesting materials. 
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